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It is shown that for a generic electrovacuum spacetime, electromagnetic radiation produces vor- 
ticity of worldlines of observers in a Bondi-Sachs frame. Such an effect (and the ensuing gyroscope 
precession with respect to the lattice) which is a reminiscence of generation of vorticity by grav- 
itational radiation, may be linked to the nonvanishing of components of the Poynting and the 
super-Poynting vectors on the planes othogonal to the vorticity vector. The possible observational 
relevance of such an effect is commented. 
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I. INTRODUCTION 

Nowadays high-accuracy experiments in space can test 
relativistic gravity. These endeavours surely should lead 
to direct detection of gravitational waves Particu- 
larly interesting are the frame dragging observations us- 
ing cryogenic gyroscopes in Gravity Probe B searches 
[3] and from data provided by LAGEOS satellites 
These successful experiments have fundamental implica- 
tions for general relativity since frame dragging is one of 
the most conspicuous general relativistic effect 0]. Fur- 
thermore, they are also very important in astrophysics 
because the frame dragging has been invoked as a mech- 
anism to drive relativistic jets emanating from galactic 
nuclei Q. Dual jets are expected for inspiral supermas- 
sive binary black hole systems which produce prodigious 
quantities of gravitational waves and energetic electro- 
magnetic events Q , which in some way have to be closely 
related. Less impressive outcomes but of utmost impor- 
tance are the intense electromagnetic outbursts from hy- 
perenergetic phenomena such as collapsing hypermassive 
neutron stars @ and Gamma Ray Bursts Q. 

Now, it is a well established fact that gravitational ra- 
diation produces vorticity in the congruence of observers 
with respect to the compass of inertia [l0[-[IIj|- Since 
the vorticity vector decribes the proper angular veloc- 
ity of the compass of inertia (gyroscope) with respect to 
reference particles [l5[, it is clear that a frame dragging 
effect is associated with gravitational radiation. 

In [bH it was further shown that such a vorticity is 
closely related to the non-vanishing of super-Poynting 
vector components on the plane orthogonal to the vor- 
ticity vector. Later it was shown that the vorticity ap- 
pearing in stationary vacuum spacetimes is also related 



to the existence of a flow of superenergy on the plane 
orthogonal to the vorticity vector [l6| . 

The rationale to link vorticity and the super-Poynting 
vector comes from an idea put forward by Bonnor in or- 
der to explain the appearance of vorticity in the space- 
time generated by a charged magnetic dipole [13] • In- 
deed, Bonnor observes that for such a system there ex- 
ists a non-vanishing component of the Poynting vector, 
describing a flow of electromagnetic energy round in cir- 
cles [l8[. He then suggests that such a circular flow of 
energy affects inertial frames by producing vorticity of 
congruences of particles, relative to the compass of iner- 
tia. Later, this conjecture was shown to be valid for a 
general axially symmetric stationary electrovacuum met- 
ric [H. 

In [l3[ it was suggested for the first time that a similar 
mechanism might be at the origin of vorticity in the grav- 
itational case, i.e. a circular flow of gravitational energy 
would produce vorticity. However due to the well known 
problems associated to a local and invariant definition of 
gravitational energy it was unclear at that time what ex- 
pression for the "gravitational" Poynting vector should 
be used. Following a suggestion by Roy Maartens to one 
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of us (L.H.) we tried in [14[ with the super-Poynting 
vector based on the Bel-Robinson tensor |2(|-[23|]. Do- 
ing so we were able to establish the link between grav- 
itational radiation and vorticity, invoking a mechanism 
similar to that proposed by Bonnor for the charged mag- 
netic dipole. 

Our purpose with this work is to tie up an important 
loose end related to this problem, namely: Does electro- 
magnetic radiation produces vorticity?, and if so can we 
explain such vorticity through a Bonnor-like mechanism? 

We shall analyze a generic electrovacuum spacetime 
following the scheme developed by van der Burg_|2J], 
which is an extension of the Bondi formalism [25ll26l] as 
to include electromagnetic fields. After a brief review 
of van der Burg paper we shall proceed to calculate the 
vorticity, the Poynting and the super-Poynting vectors. 
Based on the analysis of the obtained results we shall 
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then answer the two questions raised above. 



The metric ([T]) can be written as follows [2£ 



II. ASYMPTOTIC EXPANSIONS FOR THE 
EINSTEIN-MAXWELL FIELD 

In [24[ van der Burg presents a generalization of the 
Bondi-Sachs formalism for the Einstein-Maxwell system. 
This has, among other things, the virtue of providing a 
clear and precise criterion for the existence of gravita- 
tional and electromagnetic radiation. Namely, if the news 
functions (gravitational and/or electromagnetic) are zero 
over a time interval, then there is no radiation (gravita- 
tional and/or electromagnetic) during that interval. 

The formalism has as its main drawback [27| the fact 
that it is based on a series expansion which could not give 
closed solutions and which raises unanswered questions 
about convergence and appropriateness of the expansion. 
However since we shall consider regions of spacctime very 
far from the source, we shall use in our calculations only 
the leading terms in the expansion of metric functions. 
Furthermore, since the source is assumed to radiate dur- 
ing a finite interval, then no problem of convergence ap- 
pears [H| . 

The general form of the metric in the Bondi-Sachs for- 
malism can be written as 1261 



ds 2 = (Vr 
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r e 



C/ 2 cosh 25 -r 2 e~ 27 jy 2 x 



cosh 28 - 2r 2 UW smh28)du 2 + 2e 2f3 dudr + 2r 2 x 
{e 2l U cosh 26 + W sinh 26)dud8 + 2r 2 {e~ 2l W x 
cosh 26 + U sinh 28) sin 9dud(f> - r 2 (e 27 cosh 28d9 2 



- 2l cosh 28 sin 2 



2 sinh 26 sin 0d0d0), (6) 



if we choose the gauge of Bondi, det(h A s) = det(q A B), 
where qab is the unit 2-sphere metric; 7, 5, U, W are 
functions of (u,r,0,<p). In the present case the space- 
time is asymptotically flat and necessarily Minkowskian 
by construction. 

The procedure follows the script established by Bondi 
ct al. [25|. Thus, four functions are assumed to be ex- 
panded as power series in negative powers of r, and pre- 
scribed on a hypersurfacc u = uq =constant, namely (sec 
IH for details) 
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i>~ 4 + ..., (7) 



5 = dr- 1 + [H + c 2 d/2 - d 3 /6]r" 3 + Kr^ + (8) 



ds 2 



V 



2 h AB U A U B ) du 2 + 2e 2l3 dudr 



2r 2 h AB U B dudx A - r 2 h AB dx A dx l 



(1) 



where all the metric components are functions of x° = u, 
x 1 = r, x 2 = 9, x 3 = (j). /3 is the expansion of outward 
null rays, V the Newtonian like potential. If we make 
a 2 + 1 foliation doing constant r, e 2/3 V/r is the lapse, 
U A the shift and h A B the 2-surface metric of constant 



u, satisfying h AB h BC 



5 A . u is a timclikc coordinate, 



constant along outgoing radial null geodesies, while r is 
a luminosity-distance parameter. 

The electromagnetic field is characterized by a skew- 
symmetric tensor F M „ giving rise the energy-momentum 
tensor T^ u defined by 



T 



1 



9^ F^Fvs, 



(2) 



being g^ v the metric tensor given by ([T]). The Einstein- 
Maxwell field equations for the electro-vacuum spacetime 
outside the source are 



Tu 7 



0. 

0, 



(3) 
(4) 
(5) 



where i? M7 is the Ricci tensor 



Note that the units used 
are 16nG = c = 1. Our purpose here is to extract general 
information from the radiative zone, at J7 + (future null- 
infinity), where the well known asymptotic expansion in 
power of r" 1 of Bondi seems to be the most convenient 
and precise. 
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+ (2E + ec + fd)r~ 6 + 



F13 esc 6 = fr- 2 + (2F + ed - fc) 



(9) 



(10) 



where all coefficients are functions of u, 9 and <p. 

Next, from (|7])- (fTU|) and a subset of field equations 
(main equations), the following expressions are obtained 
for metric and electromagnetic variables: 



-(c 2 + d 2 )r- 2 /A + 



(11) 



U = — (eg + 2c cot 6 + d<j, esc 0) \ 



+ [2N + 3(cc0 + dd e ) + 4(c 2 + d 2 ) cot 9 

— 2{c^d — cdff,) esc 9] r~ 3 

+ - {3C + 2C cot 9 -+- esc 9 — 6(cN + dQ) 

- 4(2c 2 c e + eddg + c e d 2 ) - 8c(c 2 + d 2 ) cot 9 - 4(c 2 c^ 
+cc 4> d + 2d 2 d (j ,)csc9 - (ee-^/)}r~ 4 H , (12) 



W = -(d e + 2dcot9-c,p esc 9)r~ 2 

+ [2Q + 2(ced - cd e ) + 3(cc0 + dd^) esc 9}r~ 3 
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{3(H e + 2H cot 9-C^ esc 9) - 6{cQ - dN) 



- 4{2d 2 d g + cc d + c 2 de) - 8d(c 2 + d 2 ) cot 9 + A(c<pd 2 
+cdd rjj + 2c 2 c<p) esc 9 + (fie + e/)r" 4 } + ■■•, (13) 
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V = r - 2M - [N g + N cot 9 + Q esc 9 - - (c 2 + d 2 ) 

-(cj + dg)- 4(cc e + ddg) cot 8 - 4(c 2 + d 2 )cot 2 
-(c 2 + d\) esc 2 6 + 4(c0C? - cd^) esc 6» cot 9 



+2(c^d e - cgd^) csc9 - i(e 2 + ^ 2 )]r 1 
+{ /"(/<? + / cot 61 - e c sc 6>)}r~ 2 • • 



1 



F 02 =X + (eg - e u )r- 1 - {[E + -(ec + /d)] t 
+ i(e e + e cot 6» + / esc ^)] e }r" 2 + 



(14) 



-er 2 + (e<? + e cot# + / esc 9)r 3 + (15) 



(16) 



it follows from one of the supplemcntrary conditions (see 
IH for details) 



(c* c* + -X*X*) sin 0dQd4>, (22) 



o Jo 



clearly exhibiting the decreasing of the mass function in 
the presence of news (bar denotes complex conjugate). 

Finally, the following equations derived from the sup- 
plementary conditions (Eqs.(14)-(20) in j24|) will be 
used: 

K = ~« ~ \ x * x * + \C-iC-2c* u , (23) 



esc 6F 03 = Y + (c^csc9- /Jr" 1 -{[F+-(ed- fc)] u 



+ -(e e + ecoU 



, csc( 



, csc8}i 



esc 6F23 = — n — (fg + f cot 9 — esc 9)r 1 



(17) 



(18) 



where again, all coefficients are functions of u, 9 and </>, 
and subcript letters denote derivatives. 

At this point we have nine functions of three vari- 
ables which are undetermined on the initial hypersurfacc, 
namely M, N, Q, e, fi, c„, d u , X, Y. However using the 
remaining field equations (supplementary conditions), 
equations for the u- derivatives of M, N, Q, e, /it, in terms 
of the prescribed functions, can be obtained. Hence if 
7, S, F12, F13, M, N, Q, e, \i are prescribed on a given ini- 
tial hypersurface u = constant, the evolution of the sys- 
tem is fully determined provided the four functions, rc- 
ferrred to as news functions, c u , d u , X, Y are given for all 
it. In other words, whatever happens at the source lead- 
ing to changes in the field, it can only do so by affecting 
news functions and viceversa. In light of this comment 
the relationship between the news functions and the oc- 
currence of radition becomes clear. 

Furthermore, the mass fuction m(u) which coincides 
with the Schwarzchild mass in the static case, is defined 
by 



with 



m(u) 



1 



2-rr t"K 
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M* smOdOdcj), 



(19) 



M* = M + -i{d/d9 + cot 6)(d e + 2d cot 9 - Cj, esc 9) 
- csc9d/d<t>(cg + 2c cot 9 + d 4> csc9). (20) 
Then, introducing the intermediate complex quantities: 



N* = C M* + 2c*£_ 2 c* + £- 2 «c*) - e*X*, (24) 



e* — £_i X*, 



(25) 



AC* = 2c* 2 c* u + 2c* M* + CiN* + e* X\ (26) 



4£>* = -£_ 3 (£ 2 C* + 2c* N*) 



2e* u = -C c* - 2c* X*, 



(28) 



c + id, X* = X + iY, 



(21) 



AE* = -£_ 2 (£ ie *-2c*e*), (29) 

where 

N*=N + iQ, C*=C + iH D* = D + iK, ... 

(30) 

e*=e + ifx, e*=e + if E* = E + iF, ... (31) 
and 

£ P = -{d/d9 -p cot 9 + i esc 68/ defy. (32) 

We shall next calculate the general expressions for 
the vorticity, Poynting (electromagnetic) and super- 
Poynting vectors for our system. 
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III. VORTICITY, POYNTING AND 
SUPER-POYNTING VECTORS 

The vorticity vector is defined as usual by 

W Q = 1 rj arlLX U v U h x, Vap-yS = V~9 ^a0-y6, (33) 

9 

where rj a p 1 s = +1 for a, /3, 7, 5 in even order, — f for 
a, /?, 7, 5 in odd order and otherwise; u M is the 4- 
velocity vector for an observer at rest in the considered 
frame. The absolute value of us a is denoted by f2, i.e. 

fi = |w Q w Q | 1/2 . (34) 

The electromagnetic Poynting vector is by definition 

S a =T af3 Ufj, (35) 

whereas the super-Poynting vector based on the Bel- 
Robinson tensor, as defined in [2l|, is 

P a = y la p lS E p H^u 5 , (36) 

where and H^u, are the electric and magnetic parts 
of Weyl tensor, respectively, formed from Weyl tensor 



Ca/3-yS an d its dual C a p 7 s by contraction with the four 
velocity vector, i.e. 

E a/ 3 = C^psu'u 5 , (37) 



H a f) = Cajpsu^u 8 = -rj aie sC eS Pp u' y u p . (38) 

Let us first calculate the vorticity for the congruence 
of observers at rest in ((HJ), whose four-velocity vector is 
given by u a = A~ 1 5^, where A is given by 

A = (Vr^e 20 - r 2 e 2 ^U 2 cosh2<5 

- r 2 er 2l W 2 cosh28 - 2r 2 UW sinh25) 1/2 . (39) 

Thus, (J231) lead us to 

uj a = (w u ,w r ,wV*), (40) 

where 



lo u = — {r 2 er 2fi (WU r - UW r ) + \2r 2 sinh 28 cosh 25(U 2 e 2 ' 1 + W 2 e- 2 ' 1 ) + AUWr 2 cosh 2 25] e~ 2p lr 

2 A 2 sin 9 L J 

+ 2r 2 e- 2 ^(W 2 e- 27 - U 2 e 2 ^)5 r + e 2p [e- 2f) {U sinh 25 + e'^W cosh 25)] e 

- e 2f} [e- 2/3 (W sinh 26 + e^U cosh 25]$}, (41) 
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-{2r 2 A" 2 [((?7 2 e 27 + W 2 e~ 2 '') sinh 25 cosh 26 + UW cosh 2 2<5) 7u ) + (WV 



27 



e 2/3 sin 

-{WU U - UW U )] + A 2 [A- 2 (We- 2 ^ cosh26 + U smh26)] g ^ A 2 [A~ 2 (W sinh25 + Ue 2 ^ cosh26)] cj> }, (42) 



{ A 2 e~ 20 [r 2 A~ 2 (U sinh 2<5 + We~ 2j cosh 2<5)] r - e 2f3 A~ 2 [e~ 2fi r 2 (U sinh 2(5 + e~ 2l W cosh 25)], 



2r 2 sinfl 
e ^A- 2 (e- 2 ^A%} 



(43) 



and 



1 



2r 2 sin 9 

+ A- 2 e 2 ^{A 2 e- 2 ' i )e} 



{A 2 e- 2 P[r 2 A- 2 (Wsinh25 + Ue 2 "< cosh25)] r - e 2fi A- 2 [r 2 e~ 2fi {W sinh2<5 + Ue 2 * 1 cosh2<5)], 



(44) 
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Although algebraic manipulation by hand is feasible for 
the Bondi metric, for the Bondi-Sachs one is quite cum- 
bersome. Therefore we write a Maple 15 script (avail- 
able upon request) which uses intrinsic procedures to 
deal with tensors. We proceed in two steps. First, we 
calculate and save (|37|) . ([38| and ([36| . using the general 
form ((TJ), that is, without using the metric function ex- 
pansions (fTT j) -(fT4 | . Second, we expand all the relevant 
objects separately (metric, Weyl, electric and magnetic 
parts) up to the leading order. After these two simple 
steps we were able to write the output for the super- 
Poynting. 

We use the shift vector U A = (U, W/ sin 6*) and the 
2-surfaces metric of constant u 



l-AB 



e 27 cosh 25 
sinh2<5 sm( 



sinh 25 sin 9 

e- 2 ~f cosh 25 sin 2 9 



(45) 



In our calculations we keep these auxiliary variables, U A 
and hAB-, as far as was possible. 

Thus, for the absolute value of we get 



fl = flgr 



flgSM^' 



(46) 



where subscripts G, Q&M. and EM. stand for gravita- 
tional, gravito-electromagnetic and electromagnetic. At 
this point a remark on the meaning of this notation 
is in order: "gravitational" terms refer to those terms 
containing exclusively functions M, N,Q,c u ,d u and their 
derivatives. "Electromagnetic" terms are those contain- 
ing exclusively functions e, /i, X, Y and their derivatives, 
whereas "gravitomagnctic" terms refer to those contain- 
ing functions of cither kind and/or combination of both. 
It should be clear that all functions are related through 
field equations, and therefore the established division is 
rather formal, however the splitting as indicated is useful 
for the discussion below. 

In order to get some insight let us first consider the 
axially and reflection symmetric case, i.e., d/d<f> = d = 
/ = H = K = F = Y = Q = 0. Thus, we obtain 

flgsM = -j[4c u Mc e + 2N ee + 2c u N + 3Cg u + 2N 

-4(Mc) e - (ee) u + 2(4M - c)N u 
-4cc0C u + 8M g M - 2(1 + 2c u )N cot 2 9 
+(8cM - 8M 2 - 2N e - 3c 2 )c e „ 
+2(3c u c 2 - 2c u N g - 8c u M 2 + Ac u cM 
-Nc Bu - 4cM + 3c u + Ng) cot 9], (47) 

where we observe the contribution \{ee) u , which is 
purely electromagnetic. This latter term together with 
terms proportional to N u and C u may in turn be ex- 
pressed through electromagnetic news (X*) by means of 
d21J|-((2Bl), and pg|l. Doing so we clearly identify electro- 
magnetic radiation as a vorticity source. 

For the angular contravariant super-Poynting com- 
ponents (in the axially symmetric case) we identify 
many QEM. contributions. One of them is explicit, e.g. 
Cuu(ee) u in P 9 (the complete term is too long to display 



here; see the appendix) whereas other terms appear by 
replacing u derivatives of metric variables by means of 
(|32l) . It is remarkable the following simplification 



P, 



gsM 



esc 9c uu {]xe) u . 



(48) 



Let us now get back to the most general three- 
dimensional case, where except for the leading term, 

fig = ~ [(C8u + 2c u cot 9 + d<t, u esc 9) 2 

+ (d eu + 2d u cot 9 - tyu csc 9) 2 ] 1/2 , (49) 

the next leading gravito-electromagnetic term is too long 
to write here or anywhere (the expression is available 
upon request). 

However, the important point to stress here is that 
flgSM contains electromagnetic news functions (X*) by 
means of, for instance of terms 

3sm9cj u d lt>u ee u , 3ce M ^„//z u , (50) 

which are typical and selected contributing terms to 
flgsM- 

Next, we write the electromagnetic Poynting vector as 



qu _ nu —4 , CM —5 



ar _ or —1 , or —3 
O — '-'EM 1 ' KA< 



} g£M' 



08 _ o8 -4 , 08 -5 



(51) 



(52) 



(53) 



S* = St M r- A + St M r- 5 + --- , (54) 



where 



S U £M = ^+t% 



(55) 



s g£M = (A*/ _ ee ) cot 61 _ ( e U + P- e 4>) csc 

+tife + MS% M -ee e , (56) 



S r £M =X 2 +Y 2 , 



(57) 



SgsM = - 2Xe u + MS r £M + 2Xe e - AdXY 

+2c(Y 2 - X 2 ) -2Yf u + 2Y C(P csc9, (58) 
S° £M =eX + pY, (59) 

Sg £M = -VH + X U ~ V**) c s c - (Xe - Yf) cot 9 
+MS e £M - 2e{Yd + Xc) + eeg - Xe e - fxf u 
+Yf e - ee u , (60) 
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S* = csc 9(eY-fxX), 



(61) 



Other terms are too long to display. However the im- 
portant point is that, again, terms containing electro- 
magnetic news appear in Pg£j^, as for example 



S* £M = - esc 9[(Ye + Xf) cot 9 - MS% M sin 6 

-{Xej, - Y,U + ec ) cscO + (2de + f e )X 
+{e e - 2ce)Y + ef u + (eg - e„)/i] . (62) 



cot 4 9c uu ee v 



COt QCuufHu 



and 



2 esc 9 cot 6d uu ee u , esc 5 9d uu f^ u , 



(67) 



(68) 



It is worth noticing that terms explicitly containing X* 
appear in 9 and <f> component of the Poynting vector. 

Next, calculation of the super-Poynting gives the fol- 
lowing result 



P M = (O,P r ,P e ,P ), 



(63) 



where the explicit terms are too long to be written at 
this point. 

The leading terms for each super-Poynting (contravari- 



ant) component are 




S Q = s Qgr 2 + ■ 


■ ■ + s^gsM 7 " 


- 4 + 




P u = P£r~ 4 + • • • , 




qU _ qU -4 








P r = Pgr' 4 + ■■■ , 




qr qr —4 

s o — s^gsM™ 


H 






I" for 1 P» SM r ' 




qO _ qe -5 


H 






P 4> = P$r-* + --- + Pl SM r-« + ... , 


(64) 


q4> — q<P _-5 
s j — a^QSM.' 


H 






where 




S P U = s P^r- 6 + ■ 








Pg = 2\p,(d uu d u + c uu c u ) cot 8 + c uu cg u 

~~\~{ c uud(pu ^uu^u) cscO -+- d uu dg u ]^ 




s P r = s P r g r- & + ■ 








(65) 


S P° = s P*r" 7 + • 




- 9 + 
- 9 + 




Pg =2 csc 9[2(c uu d u - d uu c lt ) cot 9 + c ull de u 




S P* = s Plr- 7 + ■ 


' ' + s Pg£M r 




duuC$u (CuuCcfiu ~t~ duud(pu^) CSC . 


(66) 

1 


where 








s^G = ^{[4(d00C0 e - MgCcj,) cos9 + 2(4c<pgc^ 


— idgd^ 


0)] cos 9 + [^(dd^ - ddgg -\ 


- M 2 + M 2 + 


4e + 


d 2 ) 


+(d^ - d ee ) 2 } csc 9 + [2(2c 0e - 3dg cos 6) (2d - 


- dee) + 4(2dc0 - ic^gdg - 


dggC^) COt 9 + 


(Ml 


— a ee 



which are typical and selected contributing terms to the 

super-Poynting components Pg^j^ anc * PgSM respec- 
tively. 

It is easy to check that the stationary case satisfy well 
known results, for instance, that of Bonnor or the Kcrr- 
Newman. Now, contributions to each relevant object are 
superior for s O and s P a with respect to the general (ra- 
diative) case, but it is not true for s S a which keeps the 
same leading terms: 



(69) 
(70) 
(71) 
(72) 
(73) 
(74) 
(75) 
(76) 
(77) 



+4(dd 6 
-d 2 aa + 



d 2 - Mg)) cos6>cot6>] sin 2 9 + [Ad^cj, + (4(dd, 
cos 9] cot 9} 1/2 



96 



MS - dd. 



Mi 



+ 4) 



(78) 



(the expression for s QgEM is too long to write here 
or anywhere; we check that is manifestly gravito- 
electromagnetic) 



qU 



qU 



(79) 



q9 



ege + cgfi csc ( 



(82) 



qU 

s°g£M 



(nf — ee) cot 9 
-(e/0 + ^e^) csc6> 
+fife - eeg + MS u eM 



qr 
s°g£M 



2 2 
Cm CSC I 



q<P 



(C0f 



'9 — EgfX) CSC( 



(83) 



(80) (coefficients for s P a are too long to display anywhere). 

The above expressions illustrate the Bonnor-like mech- 
anism for any stationary clcctrovacuum solution (see [l9[ 

(81) for further discussion). 
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IV. DISCUSSION 

We have seen that electromagnetic radiation as de- 
scribed by electromagnetic news functions does produce 
vorticity. It is important to stress that this is so even 
in the case of minimum electromagnetic degrees of free- 
dom (which is one), the reflection and axially symmetric 
case, when X ^ and Y = 0, implying that the above 
mentioned effect is generic. Since vorticity unavoidably 
produces frame dragging, we have established the link 
between these two physical effects. Furthermore we have 
identify the presence of electromagnetic news both in the 
Poynting and the super-Poynting components orthogo- 
nal to the vorticity vector. Doing so we have proved that 
a Bonnor-like mechanism is at work in this case too. 
However it is important to emphasize that in the present 
situation vorticity is generated by the contributions of, 
both, the Poynting and the super-Poynting vectors, on 
the planes orthogonal to the vorticity vector. It must 
be stressed that the mechanism to explain how electro- 
magnetic and gravitational radiation produce vorticity 
invokes the concept of superenergy (its flow), which is 
one of the most important concepts in general relativity 
involving the congruence of observers [30(. 

Before proceeding further with our discussion two com- 
ments are in order: 

• We have clearly established the link between vor- 
ticity and electromagnetic radiation, which as men- 
tioned before implies a link between electromag- 
netic radiation and frame dragging. This was the 
main goal of our work. However we have not 
calculated in detail the resulting precession rate 
of a falling gyroscope under such a circumstance, 
since it is out of the scope of this manuscript. It 
goes without saying that for an explicit experiment 
setup such a calculation should be provided. 

• It should be clear that the magnitude of the men- 
tioned effect is directly related to the intensity of 
the electromagnetic emission (Y*). 

Simulations from numerical relativity could shed some 
light to figure out how to measure the effect reported 
here. First, in the study of binary black holes dynam- 



ics near electromagnetic fields and plasmas [3l[, it was 
displayed how the system imprints characteristics on the 
two induced wavebands. In the present case, as an in- 
verse problem, the electromagnetic and gravitational ra- 
diation produce precession (on test gyroscopes) which 
has to be imprinted by the waves. Second, almost in the 
same aforementioned context, it was possible to track the 
precession of compact binaries from gravitational wave 
signals [32j , locating the frame from which the (I = 2, 
\m\ = 2) modes are maximized. We suppose that in the 
same way as the "quadrupolar-aligned" frame is located, 
the "dipolar-aligned" frame could be find from electro- 
magnetic modes. This simple method can be applied to 
the ensuing gyroscope precession, as reported in here. 

The potential observational consequences of the re- 
ported effect should be seriously considered. Indeed, as 
we mentioned before, intense electromagnetic outbursts 
are expected from hyperenergetic phenomena such as 
collapsing hypermassive neutron stars and Gamma Ray 
Bursts (see [9[ , Q and references therein) . Although we 
are not able at the present to estimate the required sen- 
sitivity of the gyroscope to measure such an effect, the 
high intensity of radiation in the above mentioned sce- 
narios leaves open the question about its detection with 
present technology. In fact, the direct experimental evi- 
dence of the existence of the Lense-Thirring effect 
brings out the high degree of development achieved in 
the required technology. In the same direction point re- 
cent proposals to detect frame dragging by means of ring 
lasers 0^ 

Finally we would like to mention that frame dragging 
produced by the so called optical vortices has been re- 
cently described in the linear regime (38j . However it 
should be observed that the effect reported here stems 
from non-linear terms, as it can be seen from ([50)1 . (p7|) 
and (|68|) . bringing out the relevance of nonlinearities in 
the general relativistic description of radiation. 
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Appendix: Expression for the axially and reflection symmetric Pqsm 

P<3£M = -^[3ceC0e u c u - ^c Bu c u cgg - 4c uu e u e - 24c uu cN u - 4c tm ee u + 14cc eu c u - 60c eu c uu M 2 - 32Mcg u c u 

— 18c uu ccg — 12ce u c U uNe — 40cg u c U uC 2 — 24Mgcc uu + 8ccg u M u — 17cce u ceeu + 12Mcg u cgg u — !6c uu Mcg 
+36Mgc uu M + 30Mcg u cc uu - 18c U uMcec u - 112c u cgcc U u - 2Ac uu Nc u + 8c uu Cg u + l2Mgc u - 8cg u Ne u 
-AcouMge - l&CuuN - \2Mc0u - 6c 2 u cg + 8c uu Ngg - 8ce u c U u - lQcgc 2 0u - 6Mecgg u + (16c u cM u - 2Ac uu c 2 
-24Mc u - 16c u C uu - A&c uu Mc u c - 120c uu M 2 c u - 2c 6u M e + 8c 6u N u - 2bcc 2 6u - 12c uu Nc 9u - 8cc 2 u 

+ 12Mc 2 eu - 8c u Mgg - 32c uu cM + 2Acgg u Mc u - l6c u Ng u - 8c 2 u Cgg + 8c uu Ng - Mc 2 u M - 4lcg u Cg Cu 
+Wc uu C u - 160c uu c u c 2 - 16c u CCgg u - 24:C uu C u Ng)cot6]. 
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